Field Effect Magnetization Reversal in Ferromagnetic Semiconductor Quantum Wells 
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We predict that a novel bias- voltage assisted magnetization reversal process will occur in Mn doped 
II- VI semiconductor quantum wells or heterojunctions with carrier induced ferromagnetism. The 
effect is due to strong exchange-coupling induced subband mixing that leads to electrically tunable 
hysteresis loops. Our model calculations are based on the mean-field theory of carrier induced 
ferromagnetism in Mn-doped quantum wells and on a semi-phenomenological description of the 
host II- VI semiconductor valence bands. 



The strong dependence of transport properties in 
metallic multilayers on magnetic configuration has in- 
spired a new field of research Q, spintronics, in which 
the electron spin degree of freedom is exploited. Semi- 
conductors that become ferromagnetic when doped with 
magnetic elements, for example GaAs doped with Mn, 
are of special interest in this field in part because of their 
compatibility with conventional semiconductor technol- 
ogy and because of the larger range of possibilities for 
the control of device functionality. Ferromagnetism has 
already been observed in such systems at relatively high 
temperatures (Tc >100K) ||] and strategies for obtaining 
Tc's larger than room temperature are of a great current 
interest in this scientific community |^,^. 

Our interest here is in the magnetic properties of mod- 
ulation doped quantum wells containing Mn local mo- 
ments. We expect |^,^ these systems to be unusual fer- 
romagnets because of the reduced dimensionality of the 
itinerant hole system and the possibility of using confine- 
ment effects and doping profiles to manipulate their mag- 
netic properties. For example, the ferromagnetic critical 
temperature can be tuned by an external electric field 
as demonstrated in recent studies of (In,Mn)As field 
effect transistors ||^. In this Letter we study their hys- 
teresis properties which will be, we predict, extremely 
sensitive to external bias voltages. The interplay between 
local moment - itinerant hole exchange coupling, Zeeman 
coupling of the localized moments to external magnetic 
field, and the coupling of an external electric field to the 
orbital degrees of freedom of itinerant holes, makes it pos- 
sible to vary the coercive field by over an order of mag- 
nitude with rather modest external electric fields. As a 
result, we predict that the magnetization orientation in 
quantum wells can be manipulated electrically without 
changing the magnetic field. Our calculations also sug- 
gest that capacitance measurements can be used to probe 
the magnetic state of biased ferromagnetic semiconduc- 
tor quantum wells. 

The system we consider is a (Cd,Mn)Te quantum well 
grown in the (001) crystal direction. Unlike III-V host 
materials, where Mn acts as an acceptor, the Mn con- 



centration and the itinerant hole density can be varied 
independently in II- VI hosts. The quantum effects dis- 
cussed here are more pronounced when only one or two 
subbands are occupied, a condition that can easily be 
achieved with sizeable Mn concentrations in II- VI diluted 
magnetic semiconductor (DMS) quantum wells. 

The Hamiltonian of a DMS system is written using an 
envelope function description of valence band electrons 
and a spin representation of their kinetic-exchange inter- 
action [Bl with d-electrons M on the Mn^+ ions: 



TC ~ Ti n 



Jpd ^ Si 



Sj S{fi ~ Ri), 



(1) 



where i labels a valence band hole and / labels a magnetic 
ion. In Eq. (^, Hm describes the coupling of magnetic 
ions with total spin quantum number S' = 5/2 to an ex- 
ternal field. Si is a localized spin, Si is a hole spin, and 
Hb is a four-band envelope-function Hamiltonian for 
the valence bands. The four-band Kohn-Luttinger model 
describes only the total angular momentum j — 3/2 
bands, and is adequate at low hole densities when spin- 
orbit coupling is large. The degeneracy between heavy- 
hole {\jz\ = 3/2) and light-hole {\jz\ = 1/2) bands at 
the F-point in the bulk is lifted by size quantization ef- 
fects in a quasi-two-dimensional system. The resulting 
heavy-light gap is larger than the Fermi energy, in the 
relevant range of hole densities, allowing only the two 
heavy-hole bands to be occupied. The heavy holes have 
their spin aligned along z-axis ((001) crystal direction) so 
that the band electron spin matrix elements get smaller 
when the magnetization tilts away from the growth di- 
rection. This leads, as discussed below, to very strong 
magnetic anisotropy with easy axes along and opposite to 
the growth direction. This anisotropy is reduced in mag- 
nitude but is still large when mixing between heavy- and 
light-hole bands is fully accounted for, as we show below. 
The exchange interaction between valence band holes and 
localized moments is believed to be antiferromagnetic , 
i.e. Jpd > 0. For (Cd,Mn)Te, Jpd « 0.06 eV nm^. 

Our mean-field theory is derived in the spin-density- 
functional framework and leads to a set of physically 
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transparent coupled equations |Tl|]. The effective mag- 
netic field seen by localized magnetic ions is the sum 
of the external magnetic field and the kinetic exchange 
coupling mean-field contribution from spin-polarized car- 
riers: 



Heff{Ri) = Hext - Jpd{s{Ri))/gy.B 



(2) 



where Hext is the external magnetic field, {s{Ri)) is the 
carrier spin density at Mn sites, and 5 = 2 is the g- 
factor of the local moments. The mean polarization of a 
magnetic ion is given by 



{S){Ri) = SBs{SH,ff{Ri)/kBT)H,ff{Ri 



(3) 



where Bs{x) is the Brillouin function and Heff{Ri) is 
the unit vector along the direction of the effective mag- 
netic field defined in Eq.(||). The itinerant-hole spin- 
density is determined by solving the Schrodinger equa- 
tion for holes which experience chemical quantum con- 
finement and electrostatic external potentials, a local- 
spin-density-approximation exchange-correlation poten- 
tial, and a spin-dependent kinetic-exchange potential 
h{z). The field h{z) is non-zero only in the ferromagnetic 
state. We assume that the magnetic ions are randomly 
distributed and dense on a scale set by the free carrier 
Fermi wavevector and that any desired density profile, 
NMn{z), can be achieved during growth. This allows us 
to take a continuum limit where 



hiz) = Jpa NMniz) {S){z). 



(4) 



Results presented in the following paragraphs are based 
on a self-consistent solution of Eqs. (|^)-(^) in the zero 
temperature limit. This approach does not account for 
disorder associated with randomness in the Mn spin lo- 
cations, which could have some importance in general, or 
for direct interactions between Mn spins on neighboring 
lattice sites [|l3| . There is however a large body of liter- 
ature, mostly concerning the magnetic field dependence 
of optical properties of paramagnetic DMS systems, that 
supports the reliability of these approximations for qual- 
itative predictions. 

The magnetization reversal properties of these quan- 
tum well ferromagnets are quite unusual as we now es- 
tablish. We first discuss the dependence of the ground 
state energy on uniform magnetization orientation. As 
explained above we expect that quantum confinement 
and valence band spin-orbit coupling should conspire 
to yield a large anisotropy energy. In Fig. Q we plot 
the change in the total band energy per particle Eii{6) 
as a function of angle 6 of the magnetization measured 
from the growth direction for several 2D densities P2D- 
We assumed a (Cd,Mn)Tc/(Cd,Zn)Te quantum well of 
the width it; = 10 nm with the valence band offset 
150 meV and uniformly distributed Mn ions of density 
Nmh = 6 X 10^° cm~^. The energy increase at non-zero 



angles, i.e., the growth direction easy-axis anisotropy, is 
consistent with experiment [l4[ |, and with the expecta- 
tions based on (001) growth quantum confinement effects 
explained above. If there were no heavy-light subband 
mixing at any orientation, the magnetic condensation en- 
ergy would vanish for = 7r/2 (the operators Sx and Sy 
vanish in the heavy-hole subspace) and the anisotropy 
energy per electron would equal JSNmu/'^ ~ 45 meV, 
much larger than the values obtained in Fig. |^. This 
limit is reached, however, only when quantum well sub- 
band splittings are much larger than JSNmu and we are 
not close to this limit. Nevertheless, the magnitude of the 
anisotropy energy we obtain is still several times larger 
than in bulk ferromagnetic semiconductors []l5| . It is in- 
teresting to note that the anisotropy energy per electron 
in these ferromagnets is nearly three orders of magnitude 
larger than in cubic transition metal ferromagnets [|l6|. 
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FIG. 1. Band anisotropy energy calculated from the 
self-consistent Hartree approximation using the four-band 
Luttinger model as a function of magnetization direction from 
the growth direction. The anisotropy energy is defined as the 
energy difference compared to the growth direction orienta- 
tion. At the lowest density only one subband is occupied while 
two subbands are occupied at the two higher densities. The 
curves fit the numerical data in a uniaxial form. 

For magnetization along the easy axis, heavy-light mix- 
ing is very weak permitting the use a single spin-split 
band model with effective mass m\\ = 0.25too and out-of- 
plane mass = 0.55mo. In this (Cd,Mn)Te/(Cd,Zn)Te 
quantum well, the spin-splitting is about three times 
larger than the splitting between the two lowest majority- 
spin subbands. This enormous exchange gap between 
majority- and minority-spin states is the origin of the pe- 
culiar magnetic reversal properties in biased DMS quan- 
tum wells that we now discuss. Consider the effect of 
a field in the positive l-direction on the state where all 
spins are polarized in the negative i-direction. Because of 
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the antiferromagneticp— d coupling, spin-up holes are the 
majority carriers at zero external magnetic field. At the 
2D hole densities of interest here, only the lowest spin- 
up subband is occupied and all spin-down subbands are 
empty, i.e., the itinerant system is fully polarized. When 
the magnetic field is applied, direct Zeeman coupling to a 
local moment competes with the local mean-field kinetic- 
exchange coupling which is proportional to the itinerant- 
hole spin-density. Since the carrier density is smaller at 
the edges of the quantum well than at the center, spin 
reversal starts from the well edges. This, in turn, creates 
an exchange barrier for the majority-spins which effec- 
tively narrows the quantum well in which they reside. 
At the same time an effective double-well potential de- 
velops for the minority spins as illustrated in Fig. ||. Note 
that because of the strong kinetic-exchange potential, the 
symmetric and antisymmetric minority-spin states in the 
unbiased double- well remain nearly degenerate through- 
out the whole meta-stable part of the hysteretic magne- 
tization curve. 



(a) 



(b) 



(c) 
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FIG. 2. Schematic diagram of potentials in a DMS QW 
with external magnetic field, < |-ffeit| < \Hr\, which is 
oriented in opposition to the local spin moment. Hr, the 
reversal field, is close to the field at which the lowest minority 
spin level crosses the majority spin energy level, (a) and (b) 
show majority and minority spin potentials without a bias 
potential, while (c) and (d) show majority and minority spin 
potentials with a bias potential applied. The broken lines 
indicate the envelope wavefunctions. 

As the external magnetic field increases, the minority- 
spin energy levels are lowered and the majority-spin en- 
ergy levels are raised. When the lowest down-spin en- 
ergy level reaches the Fermi energy Ep, holes start to 
occupy the down-spin states. Our self-consistent calcula- 
tions show that once this occurs, the magnetization rever- 
sal is rapidly completed and only the uniform down-spin 
state is stable. 



This unusual process in which magnetization reversal 
and quantum confinement effects are linked suggests that 
the hysteresis curve can be modified by applying an ex- 
ternal electric field, as shown in Fig. ^j. We can param- 
eterize the electric field between the front-gate and the 
quantum well, fixing the field at the back side of the well, 
by the total two-dimensional carrier density. We have 
calculated local moment magnetization loops for several 
carrier densities where p2D = 3 x 10^^ cm^^ corresponds 
to the balanced quantum well. The results are plotted 
in Fig. ^. Two effects contribute to the change of the 
hysteresis loop with applied electric field. First, the hole- 
density profile is compressed and moves toward one edge 
of the well, permitting an abrupt localized moment rever- 
sal throughout the depleted region of the quantum well. 
Second, the increase (decrease) of the 2D hole density 
enhances (suppresses) the kinetic-exchange coupling. 
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FIG. 3. Local moment hysteresis loops in a quasi-2D DMS 
system with different carrier densities. These curves were 
evaluated in the Hartree approximation. 

The hysteresis loops shown in Fig. ^ imply that, at a 
fixed magnetic field, the magnetization can be changed 
by applying an electric field. This possibility of nar- 
rowing or broadening hysteresis loops electrically is an 
attractive feature of ferromagnetic semiconductor quan- 
tum well systems. It is somewhat reminiscent of recent 
interest in current driven magnetization reversal in 
metallic ferromagnets and could be extremely useful if 
other obstacles (most obviously the low ferromagnetic 
transition temperatures) to the use of these systems in 
non-volatile memory devices can be overcome. 

The interplay between quantum confinement and mag- 
netization reversal in ferromagnetic semiconductor quan- 
tum wells is reflected in unusual magnetocapacitance ef- 
fects. The inverse capacitance can be written as 



c- 



c: 



(5) 
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where C^^ is the inverse geometric capacitance, and 
is the inverse 2D-channel capacitance. The geometric ca- 
pacitance depends only on the distance from the quantum 
well to the gate and on the barrier dielectric constant e, 
and does not reflect the electronic structure of the 2D 
hole gas. The inverse channel capacitance 



c: 



dp 



2D 



dp- 



2D 



z 

-P2D 

e 



Ep 



is the sum of the electrostatic term related to the density 
dependent center of mass of the quasi-2D hole system 
and a thermodynamic term originating from the concen- 
tration dependent Fermi energy. This separation is phys- 
ically sensible and convenient and can be made precise 
by the arbitrary choice of a reference position outside 
the quantum well; only the total capacitance is mea- 
surable. Fig. ^ demonstrates the similarity of electric 
and magnetic responses to the electric bias field, indi- 
cating the magnetic state of this quasi-2D system can 
be detected electrically. These curves are calculated by 
starting in the partially magnetized state at the balanced 
density, p2D = 3 x lO^^cm"^; as explained previously 
the metastable state persists to the largest reversal fields 
for balanced quantum wells. When the density is var- 
ied, metastability is lost and the magnetization reversal 
is completed. The discontinuity in magnetization is ac- 
companied by a discontinuity in the electrochemical po- 
tential difference between channel and gate that appears 
as a singularity in the capacitance. 
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FIG. 4. Difference between the electric potential at a ref- 
erence point outside the quantum well and the chemical po- 
tential of the two-dimensional electron system as a function 
of carrier density, inset: local moment magnetization as a 
function of carrier density. These curves are calculated by 
starting in the metastable state of a balanced quantum well 
and varying the quantum well hole density. 

In closing we remark that this unusual magnetization 




reversal process can be operative only if the system is 
always stable against coherent magnetization reorienta- 
tion. Wc do not expect that this stability will hold for 
wider quantum wells, but believe that it does for the ge- 
ometry considered here, as we now argue. Accounting for 
the coupling to an external field along the easy axis to 
the Mn spins, the total energy per electron as a function 
of the coherent spin orientation 6 is: 



E{e) = Ei,{e) - cos{e)g^iBSH,,tNMnaHe.t)w/p 



2D 



(7) 



where C(i?e.t) = (<„ - + <J is 

the spatially averaged local moment polarization. Our 
band anisotropy energy results, summarized in Fig. 
are accurately described by a uniaxial form, Eb{9) — 
Ksvc?{9) with K = 5meV per electron for p2D — 
3 X lO^^cm^^. If K were constant and the C{Hext) 
were equal to one, the coherent rotation field would be 
He = 2K/{gp,BSNMnw/p2D) ~ 17mT. However ( is re- 
duced in the spatially inhomogeneous state and K will 
be increased by the effectively narrower quantum well 
in which the majority spin electrons are confined. These 
two effects will combine to make the unusual reversal pro- 
cess discussed in this paper an operative one for quantum 
wells narrower than ~ lOnm. 
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